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Our  investigation  within  this  project  was  focused  on  two  major  areas:  transfer  blocking  and  rejec¬ 
tion  blocking  protocols  in  queueing  networks  with  finite  capacities. 

1.  Networks  with  Transfer  Blocking 

In  this  type  of  networks  blocking  occurs  when  a  job  after  completing  its  service  at  a  station  cannot 
proceed  to  its  destination  which  is  full.  The  job  resides  in  the  server  of  the  source  station  and  waits 
there  until  a  place  becomes  available  in  the  destination  station. 

1.1.  Deadlocks 

On  page  5,  equation  (2),  of  the  proposal  [1]  we  stated  that  deadlocks  may  occur  in  this  type  of 
networks.  We  proved  this  statement  formally  in  [2],  In  [2]  we  also  gave  an  algorithm  for  the  distribu¬ 
tion  of  total  buffer  capacities  to  stations  such  that  no  deadlock  will  occur.  Additionally,  we  introduce  an 
algorithm  which  automatically  finds  cycles  in  so-called  cacti  networks.  This  paper  was  submitted  to 
" QUESTA :  Queueing  Systems:  Theory  and  Applications "  in  October  1987  for  publication.  The  paper  is 
accepted  in  June  1988  and  will  be  published  soon. 
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1.2.  Throughput  Optimization 

As  we  already  observed  in  the  investigation  of  queueing  networks  with  blocking  that  the 
throughput  is  a  non-decreasing  function  of  the  number  of  jobs  [3],  i.e.,  the  blocking  events  have  the 
effect  of  violating  the  throughput  results  which  were  mentioned  on  pages  7,8  and  9  of  the  proposal  [1]. 
Two  questions  arose  from  this  observation: 

i)  How  to  distribute  the  total  buffer  capacity  to  the  stations  such  that  no  deadlock  will  occur 
and  a  maximum  (optimal)  throughput  will  be  achieved? 

ii)  Given  the  buffer  capacity  of  each  station  in  the  network.  How  to  select  the  total  number  of 
jobs  in  the  network  such  that  the  throughput  will  be  maximum  (optimal)? 

To  answer  these  questions  first  we  assumed  that  all  stations  have  infinite  capacity  and  derived  new 
formulas  for  optimal  throughput  and  response  times  based  on  the  well-known  mean  value  analysis 
approach  [4].  Then  in  [3]  we  found  necessary  and  sufficient  conditions  for  buffer  allocation  in  cyclic 
networks  with  blocking  such  that  an  optimal  throughput  will  be  achieved. 

1.3.  Norton’s  Theorem  Application 

Our  major  breakthrough  for  this  proposal  [1]  was  that  we  found  a  very  efficient  solution  for  the 
proposed  problems  stated  on  pages  5  through  11.  The  entire  solution  is  described  in  detail  in  [6].  The 
execution  of  the  application  of  Norton’s  theorem  on  blocking  queueing  networks  as  stated  on  pages  6 
and  7  of  [1],  needs  4  steps  for  the  analysis,  (see  page  7  of  [6]): 

i)  Construction  of  the  subnetwork  T 

ii)  Throughput  Analysis  of  the  Subnetwork 

iii)  Construction  of  the  Phases  for  the  Selected  Station 

iv)  Analysis  of  the  Two-Station  Network 

As  mentioned  in  the  previous  report  for  (AFOSR  -87-0160),  the  throughput  algorithm  suggested  in 
the  proposal  [1]  (on  pages  7,8,9)  and  used  in  step  ii)  above,  has  extensively  been  studied  for  its  accu¬ 
racy.  The  algorithm  is  described  in  detail  in  [3]. 

Instead  of  determing  the  capacity  of  the  composite  station  as  proposed  on  page  9  in  [1],  we  solved 
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that  problem  in  a  different  way.  We  replace  the  subnetwork  by  a  composite  station  with  infinite  capa¬ 
city.  However,  the  blocking  events  at  selected  station  are  neglected  due  to  some  full  stations  in  the  sub¬ 
network.  Therefore,  we  modify  the  service  mechanism  of  the  selected  station  such  that  all  delays  a  job 
might  undergo  due  to  blocking  events  could  be  represented.  By  this  modification  we  obtain  a  multi¬ 
phase  server  representing  all  blocking  delays  in  the  selected  station.  We  intorduce  an  iterative  technique 
to  determine  the  parameters  such  as  the  branching  probabilites  and  delay  service  times  for  the  multi¬ 
phase  server  of  the  selected  station.  Finally,  the  entire  network  was  reduced  to  two-station  network  as  it 
was  planned  in  the  proposal  [1].  The  reduced  two-station  network  containing  the  selected  finite  capa¬ 
city  station  with  multiphase-server  representing  the  blocking  delays  and  the  composite  (flow-equivalent) 
infinite  capacity  station  representing  all  other  stations  in  the  originally  given  network  is  analyzed  by  a 
numerical  technique  as  discussed  in  [6], 

The  results  of  our  technique  have  been  validated  by  executing  several  examples  and  comparing 
them  with  simulation  counterparts  which  are  obtained  by  IBM-RESQ  package.  The  validation  study 
shows  very  good  accuracy  of  our  results. 

2.  Extended  Parametric  Analysis  of  Queueing  Networks  with  Blocking  (PART  II) 

In  the  first  step  we,  in  collaboration  with  colleagues  from  the  University  of  Erlangen -Nuernberg, 
applied  the  extended  parametric  analysis  concept  on  queueing  networks  with  infinite  capacities  and 
implemented  a  parallel  processing  on  a  multiprocessor  system  with  pyramid  type  architecture  [7].  The 
experiments  showed  that  parallel  cumputing  provides  some  advantages  regarding  storage  space.  The 
synchronization  between  processors  take  some  time  which  reflects  in  run  time  of  the  programs.  An 
open  research  question  is  how  to  optimize  the  run  time  as  well  as  how  to  optimally  decompose  the  net¬ 
work  for  parallel  computation. 

As  proposed  in  section  4  in  [1],  we  also  applied  the  extended  parametric  analysis  concept  on 
queueing  networks  with  blocking  where  the  blocking  network  can  arbitrarily  be  partitioned  into  disjoint 
subnetworks.  A  modified  and  extended  version  of  the  algorithm  presented  in  [6]  can  be  applied  solving 
the  problem  of  parallel  analysis  of  the  subnetworks.  We  briefly  outline  the  major  steps  of  the  method 
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using  the  example  from  section  4  in  [1].  The  network  to  be  analyzed  is  given  in  Figure  1. 


Figure  1. 

We  decompose  the  network  into  four  subnetworks  SNW /. 

SNW  i  =  {1,2};  SNW2  =  {12}  ;  SNW 3  =  {3 , 4 , 5 , 6 , 7  ,  8 , 9  }  ;  SNW  x  =  {10  ,  11  } 

Each  subnetwork  is  analyzed  by  shortening  all  stations  in  other  subnetworks,  i.e.,  the  service 
times  of  the  stations  not  belonging  to  a  subnetwork  are  set  equal  to  zero.  The  analysis  of  the  subnet¬ 
works  is  done  as  described  in  [6]  and  load-dependent  throughput  values  Xj(k )  are  obtained  for  each 
subnetwork  SNWj  (for  j  =  1,2,3, 4).  As  pointed  out  in  the  proposal  the  analysis  of  each  subnetwork  is 
independent  from  other  subnetworks.  Therefore,  this  part  of  the  algorithm  can  be  carried  out  in  parallel 
and  speeding  up  the  total  computation. 

In  the  next  step  we  construct  a  composite  station  for  each  subnetwork  SNWj  with  load-dependent 
service  rates  p.r;(k)  which  are  obtained  by  setting  \iCj{k)  equal  to  the  throughput  values  Xj(k).  The 
composite  stations  arc  composed  to  a  closed  queueing  network  according  to  the  configuration  of  the 
subnetworks  in  the  originally  given  network.  For  the  example  given  in  Figure  1  we  obtain  a  serially 
switched  network  as  given  in  Figure  2. 


comp,  center  4  comp,  center  3 


Figure  2. 

The  blocking  events  which  occur  during  the  transition  of  jobs  traversing  from  one  subnetwork  to 
another  are  considered  by  adding  delay  phases  to  the  load-dependent  service  stations  in  Figure  2.  The 
construction  of  the  delay  phases  is  done  as  described  in  [6].  The  application  of  our  method  to  this 
example  produces  a  network  given  in  Figure  3.  Note  that  the  phase  construction  in  Figure  3  is 
simplified  in  order  to  make  a  graphical  representation  possible. 


comp,  canter  1  comp,  center  4 


comp,  center  2  comp,  center  3 

Figure  3. 

The  final  step  of  the  algorithm  consists  of  an  iteration  between  the  network  of  Figure  3  and  the 
analysis  of  each  station  of  the  originally  given  network  in  isolation  (as  a  simple  M IM II  station  with 
finite  capacity).  Trying  to  apply  a  similar  method  as  described  in  [6]  to  this  type  of  network  we  had  to 
solve  the  following  problems: 

The  multi-phase  servers  of  the  stations  in  Figure  3  contain  load -depen  dent  phases.  We  did  not 
encounter  this  situation  in  [6]  since  the  phase  construction  was  only  done  for  the  selected  station  and 
not  for  a  composite  station  representing  a  subnetwork.  This  situation  is  solved  as  follows.  As  in  [6]  for 
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the  selected  station  wc  make  the  multi-phase  servers  computationally  tractable  by  reducing  the  service 
and  delay  phases  of  each  station  to  a  Coxian  server  with  two  phases.  The  load-dependence  is  considered 
by  a  computation  of  the  Coxian  server  under  all  possible  loads.  For  the  given  example  the  following 
calculations  are  necessary: 

fiji(t)  =  2  -flj(Jt) 

M*>  =  */(*>•  Ay  (*> 


*/!<*)  = 


1 

2  •«/(*) 


for  j  =  1,  2,  3,  4  ;  k  =  1,2,...,* 

where  K  is  the  total  number  of  jobs  in  the  network.  The  values  for  p.;(&)  and  £j(k)  are  obtained  by 
calculating  mean  value  and  coefficient  of  variation  of  the  multi-phase  servers.  This  is  straightforward 
since  the  phases  are  exponentially  distributed  and  independent  from  each  other.  As  a  result  we  obtain  a 
load-dependent  server  with  a  Coxian-2  service  time  distribution  for  each  composite  station.  From  [8] 
we  know  that  this  representation  of  a  general  service  time  distribution  provides  good  results  as  long  as 
the  squared  coefficient  of  variation  for  each  cf(k)  are  greater  than  0.5.  In  [6]  only  one  subnetwork  was 
constructed  for  the  analysis.  That  one  subnetwork  was  then  analyzed  by  the  numerical  method.  In  this 
case,  however,  the  numerical  method  is  not  applicable  because  of  the  dimension  of  the  state  space.  By 
applying  the  load-dependent  method  of  Marie  [8]  we  obtained  very  accurate  results  for  performance 
measures. 

The  applicability  of  the  extended  parametric  analysis  has  numerous  advantages.  We  mentioned  the 
acceleration  of  computing  time  by  carrying  out  subnetwork  computations  in  parallel.  Furthermore, 
some  existing  algorithms  for  very  large  finite  capacity  queueing  networks  with  have  overflows  during 
the  computation.  Using  the  extended  extended  parametric  analysis  large  networks  can  be  decomposed 
into  subnetworks  and  therefore  they  can  become  tractable  for  the  analysis. 

An  application  of  the  extended  parametric  analysis  is  given  in  [9].  There,  we  investigate  the  per¬ 
formance  of  internetwork  traffic  in  a  communication  network  where  Local  Area  Networks  (LANs)  are 
connected  using  a  Long  Haul  Network.  We  abstract  from  the  individual  structure  of  a  LAN  consisting 
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of  a  possibly  large  lumber  of  nodes  by  defining  each  LAN  as  a  subnetwork  which  is  -  applying  the 
described  method  •  reducable  to  one  composite  station.  Therefore,  we  are  able  to  parametrize  each 
LAN,  thus  allowing  us  to  focus  our  study  on  the  internetwork  traffic  without  neglecting  intranetwork 
traffic  in  the  LANs. 

3.  Rejection  Blocking 

Within  this  project  we  also  attacked  queueing  networks  with  rejection  blocking.  In  this  case  block¬ 
ing  occurs  when  a  job  after  completing  its  service  finds  its  destination  full.  The  job  is  rejected  and  goes 
immediately  back  to  the  server  of  the  source  station  and  receives  another  round  of  service.  This  is 
repeated  until  a  space  becomes  available  in  the  destination  station.  This  type  of  blocking  protocol 
appears  in  communication  networks.  For  example,  in  a  token  ring  network  suppose  a  source  node  has 
the  token  and  wants  to  send  a  message  to  a  destination  node.  If  the  destination  node  does  not  receive 
the  message,  the  value  of  the  token  remains  unchanged  and  the  token  comes  back  to  the  source  node 
which  then  assumes  that  its  message  has  not  been  received.  The  source  node  tries  so  many  times  until 
the  message  will  be  received.  This  situation  can  be  modeled  by  a  queueing  network  with  rejection 
blocking. 

In  [10]  we  give  an  algorithm  for  computation  of  performance  measures  in  reversible  networks. 
We  also  analyze  non -reversible  queueing  networks  in  [10]  where  the  analysis  is  based  on  the  duality 
concept.  Based  on  this  concept  we  give  exact  solutions  for  performance  measures.  This  paper  [10]  has 
been  presented  at  the  international  Workshop  "Queueing  Networks  with  Blocking"  in  May  1988  and  it  is 
appeared  in  the  Proccecdings  of  that  conference. 

In  [II]  we  studied  queueing  networks  with  mixed  exponential  and  non-exponential  parallel 
queues  with  interdependent  service  capacities  and  finite  common  pools  and/or  accessibility  constraints. 
We  showed  that  this  type  of  networks  have  exact  solution. 

Note  that  our  basic  concept  of  parametric  analysis  as  well  as  extended  parametric  analysis  can 
easily  be  applied  with  some  minor  modifications  also  on  networks  of  queues  with  rejection  blocking. 
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